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[57] ABSTRACT 

Nonlinear control method is provided for a closed-loop 
trajectory in a control system of the form x=/(x)+u 
through a sliding surface, s, chosen such that 
s=e-f-Xe=:, where e is the trajectory, and X is a positive 
constant using a control law of the form 

A, g-i, _ 

u = x d~ a-jj-eP* e-J(x) 

which exploits terminal attractors of the form 



where a is a constant greater than zero, fi n , fid—(2i-\- 1), 
where i belongs to the set of positive integers chosen for 
fin and fid, and fid> fin for convergence in finite time. 
For a system in which an initial S/is zero, a control law 
is used of the form 

r Pn _i Srt 

U = x d - a € - ys& - fpc), 

which yields a control equation 

S/i 

p = s + = o 

where 6„, ^=(2/+l), wherein i belongs to the set of 
positive integers chosen for d n and and 8d>8 n for 
retaining convergence in finite time. 

5 Claims, 3 Drawing Sheets 
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SLIDING MODE CONTROL METHOD HAVING 

TERMINAL CONVERGENCE IN FINITE TIME 

ORIGIN OF INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract, and is sub- 
ject to the provisions of Public Law 96-517 (35 USC 
202) in which the contractor has elected not to retain 
title. 

TECHNICAL FIELD 

The invention relates to robotic trajectory control in 
the presence of high frequency unmodeled dynamics 
enhancing the control performance of sliding control 
described by H. Asada and J. E. Slotine, Robot Analysis 
and Control , John Wiley and Sons, New York, N.Y., 
1986, pp. 139-162, by introducing terminal conver- 
gence, thus developing control laws based upon a new 
class of sliding modes referred to herein as terminal 
sliders. 

BACKGROUND ART 

Many future robot systems may be required to oper- 
ate in environments that are highly unstructured with 
varying dynamical properties, and are active, i.e., pos- 
sessing means of self-actuation. The development of 
control technologies for unpredictable environments is 
a critical first step in efforts to design autonomous ro- 
botic systems. Specifically, algorithms devised for such 
a purpose must exhibit (i) robustness to parametric un- 
certainties in dynamic models and (ii) the ability to 
adapt rapidly enough to parametric variations in order 
to insure operational performance. Although a signifi- 
cant volume of literature may be found on the problems 
of robust control and adaptive control, many issues 
pertinent to robust and adaptive control algorithms for 
large magnitude and high frequency parametric uncer- 
tainties remain unresolved. 

While positive displacement (PD) control schemes 
with gravity compensation have been shown to be glob- 
ally stable for setpoint control, they cannot guarantee 
stability in the presence of parametric uncertainties. 

Proportion, integral and differential (PID) control 
with gravity compensation is globally stable for setpoint 
control, even in the presence of uncertainty, but the 
resulting closed loop system is rather sluggish. This 
could potentially generate large transient errors, caus- 
ing tool and robot damage during contact operations. 
Compliance and stiffness control schemes guarantee 
local stability and function adequately in quasi-static 
contact situations with well characterized stationary 
environments. 

Impedance control has been applied to a wide spec- 
trum of contact applications. [N. Hogan, “Impedance 
control: An approach to manipulation: Part I — Theory, 
Part II — Implementation, Part III — Applications,” 
Journal of Dynamic Systems, Measurement and Control , 
1985, pp. 1-24] The implementations are, however, 
limited to static environments. A detailed stability anal- 
ysis of impedance control is available in N. Hogan, “On 
the stability of manipulators performing contact tasks,” 
IEEE Journal of Robotics and Automation , Vol. 4, No. 6, 
December 1988. 

Many robust control schemes have been proposed 
and implemented for pure motion and compliant motion 
control of robots and have demonstrated the ability to 
stabilize closed loop behavior in the presence of 


2 

bounded parametric uncertainties. In addition, recent 
progress in adaptive control has enabled on-line accom- 
modation of unknown robot mass and inertial parame- 
ters and payloads. [J. J. Slotine, W. Li, “Adaptive ma- 
5 nipulator control: A case study,” IEEE Transactions on 
Automatic Control Vol. 33, No. 11, November 1988, pp. 
995-1003] 

There are a number of applications in which the fore- 
going state-of-the-art control technologies would fall 
10 short of prescribed system performance. Consider two 
scenarios relevant to space applications: (1) robot- 
assisted extravehicular activity and (2) autonomous 
fresh sample acquisition during precursor science [S. T. 
Venkataraman, N. Marzwell, “Telerobots for robust 
15 space operations,” (to appear) Proceedings SPACE 92, 
Denver, Colo., 1992]. 

In the first scenario, voice-activated robots would 
assist astronauts by fetching and returning tools, helping 
transport large and heavy payloads, etc. In the second 
20 scenario, robotic elements interact with soil and rock 
media to extract core samples for scientific analysis. 
Note, however, that in both scenarios, the environment 
dynamics is partially or completely unknown. In addi- 
tion, since astronauts are active, i.e., capable of self- 
25 actuation, and environments are unpredictable, e.g., 
rocks could contain unknown crevices, shear planes and 
cracks, the rate at which the robot control law accom- 
modates the uncertainties must be explicitly controlla- 
ble. Thus, the robot control law needs to be adjusted 
30 according to the rate of change in the environment. The 
following table summarizes the environmental effects in 
the two problem areas. 


Application 

Unmodeled Phenomena 

Effect on Dynamics 

Robotic 

Cracks in Rocks 

Unstructured Dynamics 

Coring 

Holes/V esicles in Rocks 

High Frequency 


Pressure Variation in 

Variations 


Regolith 

High Amplitude 
Variations 

Astronaut 

Unpredictable Human 

Unstructured Muscular 

Assistant 

Movements 

Dynamics 


Varying Muscle Stiffness 

Unpredictable Muscular 
Actuation 


45 Further, in both of the above applications, the actual 
dynamic characteristics of the environments are highly 
nonlinear. For example, research in biomechanical sys- 
tems [J. M. Winters, L. Stark, A-H. Seif-Naraghi, “An 
analysis of the sources of musculoskeletal system impe- 
50 dance,” Journal of Biomechanics, Vol. 21, No. 12, 1988, 
pp. 1011-1025; J. M. Winters, L. Stark, “Estimated 
mechanical properties of synergistic muscles involved 
in movements of a variety of human joints,” Journal of 
Biomechanics , Vol. 21, No. 12, 1988, pp. 1027-1041], has 
55 indicated that the linear spring, dashpot model used in 
equation (8) of A. V. Hill, “The heat of shortening and 
the dynamic constants of muscle,” Proceedings Royal 
Society , Vol. 126B, pp. 136-195, would be inadequate. 

Models consisting of exponential and higher order 
60 polynomial terms have been suggested for the viscoelas- 
tic properties of the muscles, and their contractile dy- 
namics have been modeled using first-order ordinary 
differential equation (Winters and Stark, supra]. Fur- 
ther, the strong interdependences between a human 
65 neuronal impedance (input impedance of the human 
muscular structure), his mechanical impedance [Win- 
ters, Stark and Seif-Naraghi, supra], coupled with the 
variations caused by fatigue [M. Pousson, J. Van Ho- 
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ecke, F. Goubel, “Changes in elastic characteristics of 
human muscle induced by eccentric exercise,” Journal 
of Biomechanics , Vol. 23, No. 4, 1990, pp. 343-348] 
cause the environmental models to be nonautonomous. 

Similarly, environmental characteristics of rocks and 5 
regolith are extremely complex. For example, the mod- 
els suggested in D. S. Rowley, F. C. Appl, “Analysis of 
Surface Set Diamond Bit Performance,” Society of Pe- 
troleum Engineers Journal, September 1969, pp. 
301-310, for diamond coring suggests a nonlinear de- 10 
pendence of the normal thrust force on the rock hard- 
ness, characteristics of the diamond matrix, drill diame- 
ters and the drilling rate. In addition, the relationship is 
nonautonomous due to the effect of temperature on drill 
characteristics, diamond wear and chip removal. 15 

State-of-the-art PD control cannot guarantee system 
stability during the above mentioned robotic tasks. The 
sluggishness in system response with PID control could 
result in tool damage during sampling tasks and cause 
excessive human fatigue during astronaut assistance. 20 
The effectiveness of conventional compliance and stiff- 
ness control methods during autonomous sampling and 
astronaut assistance operations would be extremely 
limited given the complex nonautonomous nature of 
environment dynamics. 25 

A primary concern with sliding mode robust control 
approaches is the large switching gains required and the 
constant chattering around the sliding surface. During 
autonomous sampling, this typically causes excessive 
tool wear, sample degradation and actuator saturation. 30 
Control switching, even as a phenomenon, cannot be 
recommended for man-machine systems. Although 
chattering can be potentially eliminated by the use of 
interpolation manifolds [Asada, et al., supra, pp. 
139-157], it is not recommended for applications re- 35 
ferred to in the table above since their size cannot be 
determined accurately apriori. 

The present invention focuses on two key control 
requirements, control convergence and control robust- 
ness. When perfect model information is available, the 40 
closed loop convergence must be controllable depend- 
ing, of course, upon environmental characteristics. For 
example, during astronaut assistance, robotic tasks must 
be completed in some specified time interval. The latter 
property has been referred to as finite time control 45 
systems in the control literature. 

The second issue pertains to the development of ro- 
bust control laws that do not require high frequency 
control switches. With this motivation, a theoretical 
framework that allows terminal control convergence is 50 
developed, wherein the convergence time is finite and 
controllable. A terminal sliding mode robust control 
law is proposed to deal with model uncertainties. It is 
shown that the proposed method leads to greater guar- 
anteed precision in all control cases discussed herein. 55 

The concept of applying sliding modes to control 
emerged from earlier work on variable structure sys- 
tems, notably the work of A. F. Fillipov, “Differential 
equations with discontinuous right-hand sides,” Annals 
Mathematical Society Transactions, Vol. 42, 1964. In 60 
principle, it revolves around the choice of a control law 
that forces the closed loop system behavior to be identi- 
cal to a sliding surface. Typically, the closed loop sys- 
tem dynamics represents controlled system error behav- 
ior. As a result, one can model closed loop behavior 65 
through an appropriate choice of sliding surfaces. If a 
sliding surface is chosen such that s=e-f-Xe=0, where e 
is the trajectory error and X is a positive constant 


4 

chosen by designer, then exponential error convergence 
occurs. Consider, for example, the system 

X=/(x)+u ,(1) 

where x is the system state, / is a smooth function of x, 
and u is the control. A control law of the form u=u 0 =- 
— /(x)+Xrf— Xe will result in s— 0. If the initial condi- 
tion, s(t 0 ), is zero then the system converges exponen- 
tially. When s(t o )#0, a control switch of the form K 
Sgn(s)K>0 may be added to u^, where K is a constant, 
and Sgn(s) is 1, 0, and — 1 if s is greater than zero, equal 
to zero, or less than zero (i.e., negative), respectively, to 
force the system to converge towards the sliding sur- 
face. At s=0, the switch is deactivated, and the previ- 
ous arguments apply. It can be shown that the system 
will reach the sliding surface in finite time [Asada, et al., 
supra]. 

In the presence of parametric uncertainties where 
only /(x), an estimate of /(x), is available, control law 
takes the following form: 


« - u - KSgn(s) w 

= -f+Xd — he-K Sgn(s ) 

which yields^ s=/(x)—/(x)—K Sgn(s). By choosing 
K> || /(x)— /(x) || , Lyapunov stability and conver- 
gence towards the sliding surface ,s, can be ensured. 
[Asada, et al., supra]. The closed loop system does not 
actually stay on the sliding surface, since at s=0, s=^=0 it 
chatters in the neighborhood of the sliding surface 
[Asada, et al., supra]. Conventional sliding mode con- 
trol, therefore, guarantees exponential stability with full 
model information and asymptotic stability in the pres- 
ence of uncertainties. The design of interpolation re- 
gions is typically performed off-line using bounds on 
uncertainty and the expected system response in the 
neighborhood of the sliding surface. 

The concept of terminally sliding surfaces may be 
developed from first principles and applied to control 
synthesis for nonlinear systems. The performance thus 
obtained can be compared to those with the conven- 
tional sliding mode control. To enhance convergence 
properties of dynamical systems, the concept of termi- 
nal attractors was introduced in M. Zak, “Terminal 
attractors for content addressable memory in neural 
networks,” Physics Letters, Vol. 133, pp. 218-222, 1988. 
Since then terminally sliding surfaces have demon- 
strated considerable success in neural learning [J. 
Bahren and S. Gulati, “Self-organizing neuromorphic 
architecture for manipulator inverse kinematics,” 
NATO ASI Series, Ed. C. S. G. Lee, Vol. 44, 1990]. It 
has the basic form of a cubic parabola: 

x=-x 1 (3) 

with an equilibrium point at x=0. Integrating between 

^initial and t equilibrium* 


{^equilibrium — t initial) = " 2 * 


(4) 


This implies that Equation (4) settles into equilibrium in 
finite time. For an additional discussion of terminal 
sliders in coping with variations in unstructured envi- 
ronments, see S. T. Venkataraman and S. Gulati, “Ter- 
minal Sliding Modes: A New Approach to Nonlinear 
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Control Synthesis,” Fifth International Conference on 
Advanced Robotics, Pisa, Italy, Vol. 1, pp. 443-448, Jun. 
19-22, 1992, which by this reference is hereby made a 
part hereof. This property has also been applied for 
finite time control of distributed parameter systems. [S. 
Jayasuriya and A. R. Diaz, “Performance enhancement 
of distributed parameter systems by a class of nonlinear 
control,” Proceedings Conference on Decision and Con- 
trol, Los Angeles, Calif., December 1987, pp. 
2125-212 6]. The better convergence results from in- 
creased local stability. A detailed discussion on terminal 
attractor may be found in Zak, supra. 

The most general form for a first-order terminal at- 
tractor would be x+X(x)=0, where X is bounded for a 
bounded x Sgn(X)=Sgn x and 


dX 

dx 


Such systems are Lipschitzian in any error, e, neighbor- 
hood of the equilibrium point, but are non-Lipschitzian 
at the equilibrium point itself. An intuitive argument 
about the dynamic behavior of such systems is provided 
in M. Zak, “Cumulative effect at the soil surface due to 
shear wave propagation,” ASME Journal of Applied 
Mechanics , Vol. 50, 1983, pp. 227-229. To analyze such 
systems using Lyapunov methods, postulate the follow- 
ing: Given a dynamical system of the following form 
x+X(x)=0, and a Lyaponov function candidate V(x), 
where V is bounded for founded 
x||V(x=^0)||>0,||V(x=0)||=0,if 


V+& 0=0 


(5) 


where g is a smooth function of V, such that V has the 
terminal attractor property described above, the dy- 
namical system is terminally stable. For example, the 
system 


curve 1. Curves 3 and 4 portray the effect of attractor 
gain a. 

STATEMENT OF THE INVENTION 

5 An object of this invention is to provide robust non- 
linear controllers for robotic operations in unstructured 
environments based upon a new class of closed loop 
sliding control methods, sometimes denoted terminal 
sliders, where the new class will enforce closed- loop 
10 control convergence to equilibrium in finite time. Im- 
proved performance results from the elimination of 
high frequency control switching previously employed 
[Asada, et al., supra] for robustness to parametric uncer- 
15 tainties. Improved performance also results from the 
dependence of terminal slider stability upon the rate of 
change of uncertainties over the sliding surface rather 
than the magnitude of the uncertainty itself for robust 
control. Terminal sliding mode control also yields im- 
20 proved convergence where convergence time is finite 
and is to be controlled. A further object is to apply 
terminal sliders to robot manipulator control and bench- 
mark performance with the traditional computed 
torque control method and provide for design of con- 
25 trol parameters. 

These and other objects of the invention are achieved 
in a control system 

x-/(x) + u, 

30 for example, where x is the system state and u is the 
control for a closed loop trajectory through a sliding 
surface, s, chosen such that s— e+\e=0, where e is the 
trajectory error (x— xj), x^ is the desired trajectory, and 
X is a positive constant. The invention uses a control law 
35 of the form 


&n 


1 

jr 


u — X/j — 1 e - fix) 


40 


would be terminally stable since 
6 

implies 


K= “T y ' r - 


In this invention, terminal attractors of the form 


45 


50 


which exploits terminal attractors of the form 

0n 

x — a xfi, 

where a is a constant greater than zero, and fi n , 
fid= 1 (2i+ 1), where iel is as defined above, and fid>fin- 
In order for the initial condition s/ to be always zero, the 
following control law is used: 


a _i * n 

u = x d - e - ys& - fix) 


Pn 

i= wF* 

are exploited, where a>0 and j3„,/3<f=(2i+l), where i 
belongs to the set of selected positive integers I so that 
both fi n and fid will always be equal to odd integers, and 60 
fid> fin > For example, if the set of integers for fi n and fid 
are 1 and 2, fi n will equal 3 and fid will equal 5. The 
notation iel is used throughout with the same meaning 
and effect. 

FIG. 1 describes terminal attractors with various 65 
convergence rates. Curves 1 and 2 portray attractor 
behavior Galedon different initial conditions. Note that 
the dotted curve T is the conventional counterpart of 


55 where X is a number greater than zero selected by the 
control system designer. Substituting this equation for u 
in the system equation x=/(x)+u yields a control equa- 
tion for a surface p referred to as a dynamic terminal 
slider 

Bn 

p = s + = 0 

which convergence to a steady state s in finite time. 
These equations thus define a novel approach to a par- 
ticular class of feedback controllers which avoid the 
need for high gain control switches to redesign slider 
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7 

parameters for each trajectory in order that the initial 
condition will always be zero. 

To apply terminal sliders to control robots having 
rigid links, the dynamics of the robots may be expressed 
by 


H(x)x-+ C(x,x)x+ G(x)=t 


S d 



the guaranteed precision becomes 


where x denotes configuration variables of the robot jq 
and H(x), C(x,x) and G(x) denote the inertia, the matrix, 
Coriolis matrix and gravitational matrix, respectively. 
For robust computed torque control, 


e pt = 



(2d 


BRIEF DESCRIPTION OF THE DRAWINGS 


7= -Hix^xd+Kye+Kpe) + C(x,x)x+ G(x) 

in the equation immediately above to yield a closed- 
loop control system of the form e+KveH-K^e=0. 
When K v , K^and H are positive and definite as they are 
for robot control synthesis, the closed-loop error system 
is exponentially stable at e=e=0, which implies suc- 
cessful tracking of arbitrary trajectories with computed 
torque control using a control law of the form 

/ r -i n \ 

r = -H(x)\x d + a—- ^ c'+y(e+ae® j^J + 


which would result in a terminally stable condition. 
When model information of the system to be controlled 
is not accurate, a dynamic slider control law is used 
with a derived estimate /(x) in the form 

a & n , 5n 

u =x d - e-ys^-fix) 

When substituted for u in the control system equation 
above, the control law of the above equation yields a 
terminal slider equation of the form 


S = —ys^d + A 

where A =f(x)—f(x) . 

The equilibrium point for the terminal slider equation 
is at s=0 and 


FIG. 1 is a graph which illustrates terminal attractors 
of the form 

fin 

i 

where a>0 and fi n ,fid =(2i+ 1) , where kl, and fid>pn- 

FIG. 2 is a graph which illustrates dynamic terminal 
behavior in accordance with the present invention. 

FIG. 3 is a graph that illustrates the properties of 
terminal sliders control of the present invention. 

DETAILED DESCRIPTION OF THE 
INVENTION 

The present invention utilizes the tools discussed 
above as background to develop a new class of closed 
loop sliding mode control methods or referred to here- 
inbefore as terminal sliders for a system, such as a ro- 
botic system having a distinctive characteristic of con- 
vergence to a stable state in finite time. Their perfor- 
mance is compared with the conventional sliding mode 
summarized above. The system in Equation (1) will be 
used as a running example. As before, error e=(x— x</), 
where xj is the desired trajectory. Consider a control 
law of the form 

r Pn -i 

U = Xd - a e - J{x) 

where a,£„, and fid are defined as before. Substituting 
Equation (6) for u in the system Equation (1) , the fol- 
lowing closed loop system is derived: 

a 0n i (7) 

e + a e e — 0. 

For the above system design, a sliding surface exists 
such that 



40 


C{x,x)i + G{x) 30 


so that precision with terminal sliders becomes 


S e q — 



60 


65 


When terminal sliders are utilized with the precision of 


fin (8) 

Si = e'j -f aei&i = 0 

and the subscript i denotes initial conditions. Equation 
(8) and (7) together establish the terminal stability of the 
system in Equation (1) under the control law of Equa- 
tion (6). Under such terminal stability, the surface s is 
denoted as a terminal sliding mode surface, and control 
laws similar to Equation (6) as closed loop sliding mode 
control method or as a terminal slider control method. 
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9 

It is obvious that for feedback linearizing control with 
full model information, terminal slider control will out- 
perform its conventional counterpart. Substituting for e 
in Equation (6) in terms of e, the following control law 
is developed: 


R , W 

“~* d+a H§r e Td 

10 

For the control signal u to be bounded, for a bounded 
error e, the exponent of e must be positive. This implies 


Jn _ j_ (10) 

Pd > 2 

For the initial condition s / to be always zero, the 
slider parameters will require continuous redesign for 20 
each trajectory. In conventional sliders, this deficiency 
is overcome with a high gain switch that forces the 
system from any initial condition to the sliding surface 
[Asada, supra]. The present invention provides a new 
approach to the problem that models the behavior be- 25 
tween initial condition and the sliding surface as a dy- 
namical system. Consider the following control law 


u = xa — 



fin 

73' 


Bn 

e - ys** - Ax) 


(“> 30 


where X is a number greater than zero selected by the 
designer, and d n , 8</=(2i-bl), where fel, and dd>$n- 
Substituting Equation (11) into Equation (1) yields 


’-s-('-s-)-O-t) 

where 7)„,7)rf=(2i+ l),iel, and T]d>r]n, terminal stability 
of dynamic sliders can be easily demonstrated 


( v n I T ? it 

V S - vd - 3 , VJ 


% 3 . &n 



FIG. 2 displays the temporal behavior of dynamic ter- 
minal sliders. Curves labeled 1, 2, 3, 4, 5 and 6 are dy- 
namic sliders converging to their respective terminal 
attractors (shown by thick lines). Curves 1, 2 and 3 
show the effect of modulating the attractor gain, X, 
while curves 4, 5 and 6 illustrate the effects of manipu- 
lating the exponent 


7/1 

y d 

Curves 4, 5 and 6 have different initial conditions but 
the same starting time. 

In summary, a new approach applicable to a particu- 
lar class of feedback linearizing controllers is presented 
by Equations (6) and (11). Superior closed loop system 
performance for arbitrary initial conditions has been 
demonstrated. In addition, the need for high gain con- 
trol switches has been eliminated, thereby avoiding 
their detrimental effects during operation. 

The following applies terminal sliders to control of 
robots having rigid links. In general, the dynamics of 
the robot may be expressed as 


Bn ( 12 ) 

p — s + ys&* ~ 0 . 

40 

For any initial condition s /, Equation (12) specifies con- 
vergence to the steady state of s in finite time, after 
which the system will reach e=0 on the surface s as 
before. The surface p in Equation (12) is referred to 
herein as a dynamic terminal slider. Note that a conven- 45 
tional slider counterpart of the form s+Xs=0 would 
also provide the attraction from an arbitrary initial con- 
dition. However, the system would approach the slid- 
ing surface only exponentially. In contrast, by modeling 
the closed loop system in the form of a dynamic termi- 50 
nal slider, Equation (1 1) achieves finite time reaching of 
the slider surface from any arbitrary initial condition 
and without the application of any high gain switching. 

To demonstrate Lyapunov stability for dynamic slid- 55 
ers, consider the Lyapunov function 

2 ~7 
V=s vd , 

60 

where v n ,*u<*=(2i-{-l), id, and Vd>v„, Differentiating 
the Lyapunov function yields 

. vn Bn 

V = vd s^. 65 

vd 1 

Choosing v n and Vd such that 


H(x)x+ C&x)x+G(x)=t (13) 

where x denotes the configuration variables (joint an- 
gles, operational space coordinates, etc.) and H(x), 
C(x,x) and G(x) denote the inertia matrix, Coriolis ma- 
trix and the gravitational vector, respectively. It is well 
known that simple PD and PID control schemes [S. 
Arimoto and F. Miyazaki, “Stability and robustness of 
PID feedback control of robot manipulators of sensor 
capability,” Proc. International Symposium on Robotics 
Research , MIT Press, Cambridge Mass., 1983] provide 
global asymptotic stability for setpoint control. J. T. 
Wen and D. S. Bayard, “New class of control for ro- 
botic manipulators, Part 1, Non-Adaptive Case,” Inter - 
national Journal of Control Vol. 47, No. 5, 1988, pp. 
1361-1406, have discussed the extent of these control 
schemes. Computed torque schemes may, on the other 
hand, be applied for controlling robot trajectories [S. 
Arimoto, et al., supra]. Using these results as a baseline, 
we examine robust computed torque control, especially 
assuming to uncertainties in the gravitational and Corio- 
lis terms. 

Now consider the computed torque method for tra- 
jectory control. Given the robot dynamics in Equation 
(13), the conventional computed torque method utilizes 
the following control law: 

r- -H(x)(xd+K^+Kpe)+C(x,x)x+G(x). (14) 

Substituting this control law in Equation (13) yields a 
closed loop system of the form e+Kve+K^e=0. As- 
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suming Kp, K v and H to be positive and definite, the 
closed loop error system can be made exponentially 
stable at 6=e=0. This implies successful tracking of 
arbitrary trajectories using computed torque control 
schemes. Note that control laws that utilizeH(x^), C(x</, 5 
x</) and G(x^) have also been suggested in the past 
[Arimoto, et al., supra]. A terminal slider counterpart 
would apply a control law of the form 

/ £ Pn 8 «\ ^ 

r = -H(x) [x d + a - ^ e + y(e + ae^ J + 


C(x,x)x + G(x). 


8 n ^ W 

V = [(A - ys^) & 


-1 T 
] 





(19) 


where T denotes a matrix transpose. An important ob- 
servation that may be made is that if A is a time invariant 
constant, the closed loop system is terminally stable 
since 



driving 


o, 


This would result in 


15 


v n 


~cV vd . 


Bn 

s + y = 0 , 

20 

which is terminally stable. 

When model information is not accurate, the issue of 
control robustness must be considered. Only additive 
uncertainties are considered in the following discus- 
sions, as opposed to multiplicative uncertainties (e.g., in 25 
robot inertias). Such effects are typically reflected in the 
system’s control gains. For additional information, refer 
to Asada, et al., supra; Wen, et al., supra; and Slotine, et 
al., supra. In analyzing the example system of Equation 
(1) and extrapolating its results to control of robot ma- 30 
nipulators with uncertainties in the coriolis and gravita- 
tional terms, consider the dynamic slider control law in 
Equation (11) with a derived estimate on /(x) denoted 
by/(x) 

35 

R ^ 5n ( 16 ) 

u = Z i -a-f r <& e'-ys^-fc). 

When substituted in Equation (1), the control law in 40 
Equation (16) yields a nonhomogeneous dynamic termi- 
nal slider equation of the form 


In the preceding expression c>0, and both %and r\d are 
as defined previously. 

The stability analysis when 


may be performed as follows: It is well known that 
robots are passive systems, and that the terms I, C and 
G in Equation (13) are such that the joint velocities and 
transcendency functions involving joint positions are 
linearly separable form link lengths, masses and inertias. 
It is reasonable to assume that the model uncertainties 
arising from incorrect estimation of link lengths, link 
masses and link inertias would have a substantially 
larger effect than the uncertainties in the computation 
of transcendental functions. Consequently, A may be 
expressed using 

A=(M-M) T Y(X,X) (20) 

where Y is a matrix of joint velocities and transcenden- 
tal functions of joint positions, M is denoted a vector of 
link masses, moments and inertias, while M is the esti- 
mated value of M. Let Am denote (M— M) r . The term 


B n (17) 

i = — ys^d + A 45 

where A=/(x)— /(x). To analyze the stability proper- 
ties of the closed loop system in Equation (17), consider 
a Lyapunov function V of the form 

50 


Bn ^ 7) n (18) 

V = (A - y s 1 ^) Td 

where rfn.'qd are as defined previously, || V(A,s) || =0 
and 


Bn 

A — ys^d, 

and || V 1 1 >0 everywhere else When s is bounded, 
|| V |j is bounded since A, the uncertainty, is bounded. 
Differentiating Equation (18) results in the following: 


aA 

dt 

in Equation (19) may be expressed as 


aA 

— Aw ( 

dY 

dx a Y 

dx 

dt 

dx 

V 

ds ~ i " dx 

ds 


= Aw ( 

/ 

ay 


a Y 


s. 3 * 

aft, e + 

dx 


a Y 

dY 

dx ’■ 

* dx 


( 21 ) 


and Y are bounded and the derivatives are known to be 
continuous. Since Am in Equation (21) represents the 
effect of mass, moment and inertia it is bounded. Substi- 
tuting Equation (21) in Equation (19) the following 
equation is derived: 


-1 T 

V=[(A-y s^d) Vd ] 



Pd 

aPn 




Bn - 

i 

d d J 


(22) 
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= f(A - ***3» 


-continued 


On (26) 


J£_ jY_ ii J~7d . ^ 
as a^ *** 3x * 


Note that Equation (21) implies that as the error re- 
duces, sY/as decreases since the term dY/ax is multi- 
plied by a decaying error term. When 


. vn 

ML _ Zh. jd ■ 


in Equation (22), the closed loop system becomes termi- 
nally stable, since 




For robust control, the choice of X must be made using 
the bounds on values of various parameters, e.g., 


7 > I II Aw I! max || 


ay ay 
ax + ax 


Y. 1 1- 

‘J\\maxj—S 


20 The above condition insures terminally stable robust 
control of robots with prior knowledge on the bounds 
of uncertainty in link masses, moments, inertia, 


HI 


where c>0, and both 7) n and 7]d are as defined previ- 
ously. When 


-^Land^- 
ax u ax 


\ ds 1 


The implications of the proposed terminal slider ro- 
bust control law may be analyzed as follows. Since the 
stability properties of the system in Equation (17) de- 
30 pend on 


the system is neutrally stable and remains without 
bounded error for bounded 


35 the system is robust to configuration and/or velocity 
independent effects such as Coulomb friction. Se- 
condly, the stability of the system depends upon the 
term 


Next consider the design of parameters a, fin, fid, 
8,\n and Xd for terminally stable, closed loop system 
response using the expression 


(denoted henceforth by 


o , pn 

as ax ap n + ax 

The first step is to choose the terminal exponent param- 
eters fi n and fid- Subsequently, a choice of the terminal 
slider gain, a, may be made using: 


fid i- 


in Equation (22), rather than the magnitude of A itself, 
as in conventional sliders. This implies an accommoda- 
tion of the controller to the first order dynamics of the 
uncertainty itself. Further, as the system converges, its 
ability to accommodate 


Note that the norm of the error has not been used in 
Equation (24) since the exponent term removes the sign. 
Also, a strictly “greater than” has been used instead of 
= since 1/a must be bounded for all e. The second step 
is to choose appropriate values for 3 n and 8 d - For termi- 60 
nal stability, 


increases. This is due to the fact that the stability prop- 
erties of Equation (22) depend upon the magnitude of 


* ® n 
lZL s Td 

8 d S 




( 25 ) For values of 


A choice of y that satisfied the condition in Equation 
(25) is 


Ml < i 
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increases as || s || decreases. Knowledge of the bounds 
on 


ay 

dx 


and 


87 

ax 


can easily be obtained using kinematic information on 
degrees of freedom of the robot. Bounds on Am may be 
evaluated based upon the degree of uncertainty in the 
estimation of respective parameters. Note that the equi- 
librium point (eq) for the system in Equation (17) is at 
s=0 and 


A min 
ay 


5 This changed only by a factor of 17. Conversely, when 
terminal sliders are utilized, from Equation (28), the 
guaranteed precision becomes 


10 



(35) 


15 



where the subscript pt denotes precision of terminal 
sliders of the present invention. Equation (35) may be 
rewritten as 


Therefore, the guaranteed precision with terminal slid- 
ers becomes 


&d 

~5n 


s eq 

and 

e eq * 




(f) 


dd ££ 


20 


(28) 


(29) 


Bd Sd 


8d Bd 

~on \7Zn 




25 Typically, 

Amin 
ya 


< 1, 


30 and the exponents and 


fid _ $d 
-s — and— , 
fin o n 


( 36 ) 


Terminal slider control outperforms conventional 35 
sliders under all of the above mentioned conditions. 
Also, guaranteed precision is better with terminal slid- 
ers due to the following. In conventional sliding mode 
control [Asada, et al., supra], the interpolation region O 
is designed such that 40 


<J>-f.a4>=A 


( 30 ) 


are greater than unity. Thus, 


s d Bd 



Further, choosing 


for some a>0. This implies that at steady state, 


45 


a > 1 


fi± 

fin 


> 



(37) 


= ^ (31) 

Weq - a 

The best performance results from this system when the 5Q 
uncertainty is minimal at the closed loop systems equi- 
librium point. The best guaranteed precision is then 
given by 

Amin ( 32)55 

e PC- — 

where the subscript pc denotes precision of conven- 
tional sliders, and A min represents the minimal value of 
the uncertainty. If conventional sliders were applied in 60 
a dynamic sense (as with dynamic terminal sliders), such 
that 


p =s+\s — 0 where 

(33) 

s=e+ae=o 

(34) 


it can be insured that e p t<e pc , thereby suggesting that a 
greater guaranteed precision can always be maintained 
using terminal sliders, as compared to conventional 
sliders. 

FIG. 3 illustrated the properties of robust terminal 
slider controls. Curves A and B show the phase space 
behavior of dynamic sliders (s vs s plots). With full 
model information, from any initial state (S^i Or S51) 
the system will converge along the respective terminal 
sliders to the origin of the graphs in FIG. 3. Note that 
uncertainty term A and the term 


have also been marked along the vertical axis. Constant 
A’s are represented by the lines marked 1 and 2. When 
A is constant, the closed loop system is stable from any 
initial condition. For example, from an initial condition 
S A 0 (marked by S and A coordinates), the system will 
converge to an equilibrium S^j?. 


then epc would be equal to 



17 


5,371,669 


18 


is represented by the regions in the figure marked at L 
If the system remains in region I at all times, it will 
converge from any initial condition for all a,X>0 and 
appropriate choice of terminal exponents. For example, 
if 



is as shown by curve C, then from initial condition S 52 , 
the system will converge to some steady state Sbe- 
When 

3A 

ds 

lies in the regions marked II, the parameters a and X 
must be designed for control robustness. Let the bound 
on 



fin 5y- ! ' * V 

u = x d - CL-j^-eW e-J{x) 

5 

where a is a positive constant and fi m fid are chosen 
parameters equal to ( 2 i-f- 1 ) where i belongs to a set of 
positive integers chosen for fi n and fid so that both of 
said parameters will be odd integers and fid>fim 
10 whereby convergence to steady state is achieved in 
finite time. 

2. A closed-loop sliding mode control method for a 
system of the form x=/(x)-f u with convergence to 
equilibrium, robustness to parametric uncertainties, and 
stability upon the rate of change of uncertainties over a 
sliding surface, where x is the system state and u is the 
control for a closed-loop trajectory through a sliding 
surface s=e+Xe=0 where e is the trajectory error 
20 (x— xj), and xd is the desired trajectory, and X is a posi- 
tive constant, by applying to said system a control of the 
form 

0 fa 

25 u-Xi-a-fa-OB-'—yje-JV). 


be denoted by the line D. Note that 


where a is a positive constant and fi n , fid are chosen 
parameters equal to ( 2 i+ 1 ) where i belongs to a set of 
30 positive integers chosen for fi n and fid so that both of 
said parameters will be odd integers and fid>/3n, which 
yields a control equation 


The terminal sliders parameters must be chosen such 
that the initial value of s must be less in magnitude than 
the value Sb 0 , where S b 0 is the projection of So on the 
s axis and So is the intersection between the 


versus s hyperbolic curve and the line marked D. For 
such a design, the closed loop behavior would resemble 
the curve marked G from So to S^. Note that as the 
value of s grows smaller for a constant value of y, the 
system will be robust to increasing values of 

3A 

35 * 

i e., larger values of 

In conclusion, greater performance with sliders has 
been made possible by this invention than has been 
available with the prior art cited, which by their refer- 
ence herein are hereby made a part hereof. 

We claim: 

1. A closed-loop sliding mode control method for £ 
system of the form x=/(x)+u with convergence to 
equilibrium, robustness to parametric uncertainties, and 
stability upon the rate of change of uncertainties over a 
sliding surface, where x is the system state and u is the 
control for a closed-loop trajectory through a sliding 
surface s=e+Xe=0 where e is the trajectory error 
(x—x</), and x^is the desired trajectory, and X is a posi- 
tive constant, by applying to said system a control of the 
form 


Sn 

p — i 4- ys^ — 0 

in which an initial state s/is to be always zero and where 
X is a chosen parameter, 6 „, d d are chosen parameters 
equal to ( 2 i+ 1 ), and i belongs to a set of positive inte- 
gers chosen for and 5^ so that both of said parameters 
will be odd integers and dd>$n, whereby convergence 
to a steady state is achieved in finite time. 

3. A closed-loop sliding mode control method for a 
system of the form x=/(x)+u with convergence to 
equilibrium, robustness to parametric uncertainties, and 
stability upon the rate of change of uncertainties over a 
sliding surface, where x is the system state and u is the 
control for a closed-loop trajectory through a sliding 
surface s— e+Xe=0 where e is the trajectory error 
(x— xd), and Xd is the desired trajectory, and X is a posi- 
tive constant, when information about said systems is 
not accurate, thus requiring the use of a derived esti- 
mate /(x) for /(x), by applying to said system a control 
in the form 

o fa , fa 

u = Xd- e-ys^-fix) 

where a is a positive constant and fi n , fid are chosen 
parameters equal to ( 2 i+l) where i belongs to a set of 
positive integers chosen for fi n and fid so that both of 
said parameters will be odd integers and fid>fin, which, 
when substituted in said control method of the aforesaid 
form x=/(x)+u, yields a closed-loop sliding mode con- 
trol method of the form 


35 
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Sn 

-ys^d + A, 


20 

and fid are first chosen, and then a terminal slider gain a 
is selected such that 


in which an initial state s / is to be always zero, and 
where A =/(x)— /(x) and which yields a control equa- 
tion 


Sn 

p = i -f ys^ = 0 

where X is a chosen parameter, 8 nt 3d are as before 
chosen parameters equal to (2i+ 1) where i belongs to a 
set of positive integers chosen for d n and dd so that both 
of said parameters will be odd integers and 
whereby convergence to a steady state is achieved in 
finite time. 

4. A closed-loop sliding mode control method for a 
system as defined in claim 3 wherein said system is a 
terminally stable, closed-loop robotic system having 
rigid links and joints, said control method being carried 
out by controlling said system in such a manner that 


5 


a > 


Pd 

Pn 


1- 

e 


fin 

Id, 


after which 8 n and 8^ are chosen for terminal stability 
that satisfies a condition 


Am 


dY 


yS n 



<o 


15 where M denotes a vector of link masses, moments and 
inertias, M is the estimated value of M, and 
AAf==(M— iCl) 7 ’, where T denotes a matrix transpose, 
and Y is a motion of joint velocities and transcendental 
functions of joint positions of the robot links. 

20 5. A closed-loop sliding mode control method for a 

system as defined in claim 4, said control method being 
constrained for robust control by imposing bounds on 
values of parameters thus 


a Y 

Bs 


bY ft* ay 


25 


where Y is a matrix of joint velocities and transcenden- 
tal of joint positions, terminal exponent parameters fi n 


30 


y > 



* 
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